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ABSTRACT 

GRB 980425 and GRB 031203 are apparently two outliers with respect to the corre- 
lation between the isotropic equivalent energy Eiso emitted in the prompt radiation 
phase and the peak frequency -Epoak of the spectrum in a vFi^ representation (the 
so-called Amati relation). We discuss if these two bursts are really different from the 
others or if their location in the i?iso-E'pcak plane is the result of other effects, such as 
viewing them off-axis, or through a scattering screen, or a misinterpretation of their 
spectral properties. The latter case seems particularly interesting after GRB 060218, 
that, unlike GRB 031203 and GRB 980425, had a prompt emission detected both in 
hard and soft X-rays which lasted ^2800 seconds. This allowed to determine its i?pcak 
and total emitted energy. Although it shares with GRB 031203 the total energetics, 
it is not an outlier with respect to the Amati correlation. We then investigate if a 
hard-to-soft spectral evolution in GRB 031203 and GRB 980425, consistent with all 
the observed properties, can give rise to a time integrated spectrum with an iJpeak 
consistent with the Amati relation. 
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1 INTRODUCTION 

Amati et al. (2001) found a correlation between the energy 
emitted (assuming isotropy) during the prompt phase of 
Gamma-Ray Bursts (GRBs), -Eiso, and the frequency where 
most of this energy is emitted, Speak (the so called Amati re- 
lation). Later, Ghirlanda, Ghisellini & Lazzati (2004) found 
a much tighter correlation between -Bpoak and the collima- 
tion corrected energy E.y, for those bursts of known jet semi- 
aperture angle 6j (the so called Ghirlanda relation). These 
authors also showed that all but two bursts of known red- 
shift and -Epeak are consistent with the Amati relation. This 
was confirmed more recently by Ghirlanda et al. (2005, see 
also Bosnjak et al. 2006) with a large sample of BATSE 
bursts with pseudo-redshifts (taken from Band, Norris & 
Bonnell 2004), despite some claims of the opposite (Nakar 
& Piran 2005; Band & Preece 2005). The two outliers with 
respect to the Amati relation (and also with respect to the 
Ghirlanda relation) are GRB 980425 and GRB 031203. Both 
of them are associated with an observed supernova of type 
Ic, i.e. GRB 980425-SN 1998bw (Galama et al. 1998) and 
GRB 031203-SN 20031w (Malesani et al. 2004, Thomsen et 
al. 2004), but also GRB 030329 (i.e. Stanek et al. 2003), 



GRB 021211 (Delia Valle et al. 2002) and the very recent 
GRB 060218 (Campana et al. 2006; Plan et al. 2006) are as- 
sociated with a SN. Unlike GRB 980425 and GRB 031203, 
the other bursts so far associated with a supernova event 
obey the Amati relation. The two outliers are very close to 
Earth, with a redshift oi z = 0.0085 (GRB 980425, Tin- 
ney et al. 1998) and z = 0.106 (GRB 031203, Prochaska 
et al. 2004). Interestingly, GRB 060218 has an intermedi- 
ate redshift z = 0.033 (Mirabal & Halpern 2006; Masetti 
et. al. 2006), therefore it is a factor ~3 more distant than 
GRB 980425 and a factor ~3 closer than GRB 031203. Note 
that the other two GRBs associated with a supernova have 
z = 0.168 (GRB 030329, Greiner et al. 2003) and z = 1.01 
(GRB 021211, Vreeswijk et al., 2002). 

The aim of this paper is to discuss if these two bursts 
are real or only apparent outliers. We will explore three 
different possibilities: 

(i) these two GRBs are seen off axis, as proposed by 
Ramirez-Ruiz et al. (2005); 

(ii) what we see is the radiation surviving after having 
crossed a Thomson scattering screen; 

(iii) the real prompt emission of these bursts is subject to 
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a relatively strong hard to soft spectral evolution, and what 
we have seen is only the hardest, short-duration part of a 
much longer emission. 

Possibility i) is quite popular, and postulates that the 
two outliers are normal GRBs observed off axis, at a viewing 
angle 6^ of the order of twice the aperture angle of the jet 
9j. (Ramirez-Ruiz et al. 2005; see also other off-axis models 
by Yamazaki, lonetoku & Nakamura, 2004; Eichler & Levin- 
son 2004) The appeal of this idea is that, within the homo- 
geneous jet scenario, such outliers must exist, and should 
outnumber "normal" GRBs by a large factor, given by the 
solid angle ratios. Of course, since the received radiation for 
off-axis observers is strongly dimmed by de-beaming, the 
actual number of observable outliers depends on an accu- 
rate calculation of how both the peak energy and the total 
time integrated flux depend on the viewing angle 9^, the jet 
aperture 9j and the bulk Lorentz factor F. 

Possibility ii) postulates the existence of some scatter- 
ing material along the line of sight, following the original 
idea of Brainerd (1994) and Brainerd et al. (1998), and later 
used also by Barbiellini et al. (2004). In this case we are re- 
ceiving the transmitted flux piercing through the scattering 
screen. The received flux can be much dimmer than the in- 
cident one, and the energy decline of the (Klein-Nishina) 
scattering cross section imprints an important modification 
on the transmitted spectrum making it harder and peaking 
at a larger -Epeak- The intrinsic spectrum of the two outliers 
is then softer and brighter, and thus possibly consistent with 
the Amati relation. 

Case iii) is based on spectral evolution. For GRB 
031203, in fact, there is a strong observational evidence that 
what INTEGRAL has seen is only a part of a more com- 
plex story: besides the high energy emission above 20 keV, 
observed to a have a flat spectrum and lasting for ~30 sec- 
onds, there must be a signiflcant emission at much softer 
energies, lasting for much longer (i.e. more than a few hun- 
dred seconds, see section 4) . This has been inferred through 
the X-ray flux scattered by dust layers in our Galaxy and 
producing time variable expanding rings around the source 
(Vaughan et al. 2004). The fluence of this component may 
have exceeded the fluence above 20 keV, and this supports 
the idea that it belongs to the prompt emission. The peak 
frequency of the combined soft and hard components, al- 
beit uncertain, could well be in the range 3-10 keV, making 
GRB 031203 consistent with the Amati relation. In explor- 
ing this idea, we flnd strong similarities (and some differ- 
ences) with the behaviour of a recent bursts, observed by 
SWIFT, namely GRB 060218. Its prompt emission lasted 
for ^^3000 seconds, enabling the XRT instruments to follow 
most of it: for this burst we then have direct information 
of the high and low energy spectral components. Spectral 
fitting of the combined (XRT-I-BAT) spectra yields a peak 
frequency perfectly consistent with the Amati relation. Since 
its bolometric isotropic emitted energy is almost the same 
of GRB 031203, it is natural to consider them as twins. 

Then the question: is it possible that even GRB 980425 
had a similar behaviour? Is it possible that the prompt emis- 
sion as observed by BATSE on one hand, and by the WFC 
and GRBM instruments of BeppoSAX on the other, was 
only part of a more complex story, and that the time in- 



tegrated Epoak is much smaller than what measured in the 
first ~50 seconds? 

To determine if the two outliers are really so or they 
only appear as such, is important for several reasons. 

For instance, Paczynski & Haensel (2005) recently pro- 
posed that all "classical" GRBs are associated to the tran- 
sition from a neutron to a quark star. The two outliers, in 
this framework, should correspond to massive stars outgo- 
ing the SN Ic explosion without undergoing the transition 
to the quark star. Clearly, in this view, GRB 980425 and 
GRB 031203 are real outliers, emitting their prompt radi- 
ation with different mechanisms than the other, "classical" 
GRBs. Moreover, if these two bursts are real outliers they 
should have a flatter luminosity function with respect to 
"normal" GRBs in order to account for their lower detec- 
tion rate, as recently discussed by Liang, Zhang and Dai 
(2006). 

On the contrary, if the two GRBs are not outliers, than 
we strengthen the Amati relation, which becomes more uni- 
versal. Its physical interpretation therefore becomes even 
more compelling. 



2 HOMOGENEOUS JET SEEN OFF AXIS 

Assume that the radiation is coUimated into a cone of semi- 
aperture angle 8j and that the angle between the jet axis 
and the line of sight is Ov The fireball emits, in its comov- 
ing frame, a time integrated spectrum described by a dou- 
ble power law, smoothly joining at some frequency i/^. This 
function is a simplified version of the popular Band model 
(Band et al. 1993) composed by two powerlaws smoothly 
joined by an exponential cutoff. Calling x' = v' jv'^, the time 
integrated energy spectrum, £'{x'), is 

£'(oc') = k (1) 

where qi and Q2 represent the power laws' energy indices 
below and above v^, respectively. Since the velocity vectors 
of the fireball are radially distributed within 9j, we need 
to integrate, with the appropriate beaming factor, over the 
entire fireball surface in order to derive the energy spectrum 
observed at 9^. Calling 5 = [r(l — (3cos9)]^^ we have 

E{x,9^) = / AtP £'{x/ 5) sin 9d9 (2) 

Jmax(0, 9v-6lj) 

Within the limits of integration, the factor Atp is given by 
(see Ghisellini & Lazzati 1999) 

At/) = 2tv; 9 < 9^ - 9^ 

At/) = 7r + 2arcsin f j; 9 > 9j - 9^ (3) 

We then define -Bpcak(^v) as the maximum of the xE{x,9v) 
function. In Fig. we plot Epcak as a function of the time 
integrated fiux E{9v) = J E{x,9v)dx, having assumed 9j = 
5°, and varying 9^ between 0° and 20°. The three curves 
correspond to F = 50, 100 and 200. It can be seen that the 
^?pcak oc E^^^ behaviour (dotted line in Fig. is reached 
only for large viewing angles. In the insert, the zoom for 
small viewing angles shows that E changes by a factor 2 for 
< 6lv < ^j, while -Epcak remains approximately constant. 
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Figure 1. The peak of the observed spectrum Speak as a function 
of the time integrated flux E. Both depend on the viewing angle 
6v In the insert we show a zoom for small viewing angles, within 



(green dots) and outside 6j. We assumed 



< 6lv < 20° 



oi = 0.5, a2 = 2. The dotted line, shown for comparison, has a 
1/3 slope. The three lines have T = 50, 100 and 200, and S^^^^ = 
1.25 keV. Black points, connected with dashed lines, correspond 
to the same viewing angle 9^ = (7.5, 12.5) for the three different 
choices of F. 



The reason for the departure from the Speak oc E^^^ 
behaviour is the following (see also Eichler & Levinson 2004 
and Toma et al. 2005): for large viewing angles, there is a 
small difference between the flux received from the near and 
the far edges of the jet, which therefore contribute equally to 
the total flux. As long as this is the case, we have Speak oc 
E^^'\ As the viewing angle is decreased, the difference of 
the fluxes received from the edges increases, with the far 
edge becoming progressively negligible (with respect to the 
near edge), causing a deficit of total fiux with respect to the 
extrapolation from larger viewing angles. 

The reason of the factor 2 difference between the en- 
ergy observed at = and 6^ — 9j is the following: on axis 
observers receive most of the flux from a circle of angular 
radius l/F. At the border of the jet (i.e. for 6v ~ dj) only 
~half of the circle is inside the jet and emitting. The max- 
imum value of Speak is not observed by observers perfectly 
on-a^xis, but slightly off. To explain this effect consider the 
on-axis observer: in this case the maximum flux is again re- 
ceived by the ring of angular radius l/F, with contributions 
both from smaller angles (with larger Speak) and larger an- 
gles (with smaller Speak)- The observed Speak is found by 
integrating the spectra over the jet solid angle. For slightly 
off-ajcis observers, the contribution from large angles (on one 
side of the jet) is partially missing (lying otside the jet solid 
angle), while it appears a contributions from still larger an- 
gles from the opposite side of the jet. This latter component 
is too de-beamed to compensate for the missing one and the 
net effect is to slightly increase the observed Speak- 
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Figure 2. The "energy function" N{E) (i.e. the time integrated 
luminosity function) for the 3 cases corresponding to F =50, 100 
and 200 (upper, intermediate and lower lines) shown in Fig. [Tl 



2.1 The "energy" function 

Assume that a subset of bursts have all the same radiated 
energy Son if seen on axis. Bursts observed off-axis will be 
observed to emit an energy S, and their number is propor- 
tional to the corresponding solid angle (see Urry & Shafer 
1984 and Celotti et al. 1993 for blazars): 



iV(S,Son)dS = — = sin^ 
2n 

leading to 

\ a cos ) 



(4) 



(5) 



If bursts have an intrinsic distribution ^(Son) of energetics 
as seen on axis, than we have to integrate N{E,Eon) over 
that distribution. 

Fig. 13 shows 7V(S,Son) for the three cases shown in 
Fig. [H (i.e. F = 50, 100 and 200). The resulting function is a 
power law at small energies, and note the increase towards 
the upper end of the energy range. This corresponds to ob- 
servers within the jet aperture angle, who see almost the 
same energetics. At small energies, the slope of the power 
law is 4/3. This value has a simple explanation. For large 
viewing angles, in fact, the observed energy S oc 5^ oc 
[F(l - cose*)]-'*. Therefore (dE/dcos9)-^ oc S"^ oc S-^/^ 
which is a very good approximation of the numerical results 
shown in Fig. |2f . 



2.2 Application to GRB 980425 and GRB 031203 

If GRB 980425 and GRB 031203 were off-axis events and we 
require their on-axis energetics to satisfy the Amati relation. 



^ This results is general, as demonstrated by Urry & Shaefer 
(1984): for beaming amplification factors L oc (5P, the result- 
ing low luminosity branch of the observed luminosity function 
N{L) oc L-i-i/p. 
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then we end up to isotropic energies larger than 10^'' erg 
(even larger than suggested by Ramirez-Ruiz et al. 2004 for 
GRB 031203 using a _Eiao oc E^^^^y. approximation). 

All models invoking an intrinsic very large energy out- 
put, similar or even larger than the energy output of "clas- 
sical" GRBs, suffer from a severe problem: the two outliers 
are the closest of all GRBs and it is therefore extremely un- 
likely that such energetic bursts exist in our vicinity, for all 
reasonable luminosity functions (or energy functions). This 
issue has been discussed in Guetta et al. (2004) and with 
the detailed "trajectories" calculated in the previous sec- 
tions the problem is even worse than assuming the naive 

1/3 

iSpeak OC EJ^ bchaviour, since the steeper dependence re- 
quires larger values of Eiso to become consistent with the 
Amati relation. For this reason, we discard this possibility. 

On the other side, if jets are homogeneous, we must 
observe some of them slightly off axis. For a given sensitivity 
threshold, the amount of them depends on the specific value 
of the bulk Lorentz factor F, since larger F make the received 
flux to decrease more rapidly increasing the viewing angle. 
This can be seen in Fig. where we have marked the points 
corresponding to the same viewing angles: increasing by a 
factor 2 the bulk Lorentz factor implies to decrease by an 
order of magnitude the measured -Eiso- Fig. |5| makes this 
argument more quantitative, showing the relative number 
of expected off-axis sources with respect to on-axis ones. 

This is not meant to discard the fact that some GRBs 
can be seen off-axis, but only to discard the fact that this is 
the only reason for bringing the two outliers over the Amati 
relation. 



3 PIERCING THROUGH A SCATTERING 
SCREEN 

As already discussed, the vicinity of the two outliers argues 
against these two bursts being normal powerful bursts seen 
off-axis, (with, say, Eieo ~ 10^^ erg or more, as suggested by 
Ramirez-Ruiz et al. 2004) since any reasonable luminosity 
function would make them extremely rare in the volume cor- 
responding to their redshifts. Instead, their vicinity argues in 
favour of an intrinsically small isotropic energy and luminos- 
ity (see e.g. Guetta et al. 2004) . We then investigate wether 
these two bursts are observed through a screen of large opti- 
cal depth, which at the same time decreases their apparent 
isotropic energy and increases their apparent -Epcak. Their 
intrinsic emission (i.e. before passing through the scattering 
screen) could then be consistent with the Amati relation. 



3.1 Transmitted spectra 

The spectrum, incident on a screen of large optical depth, 
is assumed to be a smoothly joined broken power law (with 
energy indices a and (3 at low and high frequencies, respec- 
tively) as in Sec. 2, or a single power law ending with a an 
exponential cut: 



F{x) = 
F{x) = 



{x/xt)~ 



l + ix/xb^-'-' 

Fox~°' exp{—x/xc) (6) 

where x = hv / {met?) is the photon energy in units of the 
electron rest-mass energy. We assume that a < 1 and /3 > 1. 



O 
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Figure 3. Transmitted spectrum for different values of the Thom- 
son optical depth tt, as labelled. The incident spectrum has 
£^pcak = 511 keV, a = 0.5 and /3 = 2. The chosen value of a cor- 
responds to the hardest synchrotron slope produced by a cooling 
population of electrons. The dashed line, <x E^^^ corresponds to 
the hardest possible synchrotron slope produced by an electron 
distribution with a low energy cut— off, and no cooling. Harder 
spectra are difficult to reconcile with the standard synchrotron 
interpretation (Ghirlanda, Celotti & Ghisellini 2003). Note that 
we can have transmitted spectra harder than a = —0.3, even with 
incident spectra having a = 0.5, for tt > 6. 



Thus, in a xF{x) representation, the peak energy Xpeak is 
2;peak ~ Xb [(1 — a) / (f3 — l)]i/(''~") for the broken power law 
case and Xp^ak = Xc{l — ce) for the cut-off power law case. 

Along the line of sight the Thomson optical depth is 
Tt ~ aTTiAR, where n is the number density and AR is 
the typical size of the scattering region. We will consider 
values of tt even larger than one, and this is possible only 
if the scattering region is close to the burst, and becomes 
almost completely ionized by the prompt emission. Other- 
wise, we would detect very large column densities in the 
X-ray spectra of the afterglow. The Klein Nishina optical 
depth is tkn = tto"kn/o"t. The transmitted luminosity is the 
luminosity passing through the scattering material without 
sufi^ering any scattering: 



Ftix) = F(a:)exp( 



-Tkn) 



(7) 



Fig. I^shows the transmitted spectra for different values 
of Tt for an incident spectrum peaking at -Bpeak ~ 511 keV, 
with a — 0.5 and (3 = 2. Incidentally, note that the assumed 
Q corresponds to the spectrum produced by a cooling popu- 
lation of electrons, but with an initial low energy cutoff (see 
Ghisellini, Celotti & Lazzati 2000). The effect of scattering 
hardens the transmitted spectrum, to reach values harder 
than a — —1/3 (the "line of death" for the synchrotron 
interpretation, Preece et al. 1998), for tt > 6. 

The other effect imprinted on the transmitted spec- 
trum by the scattering screen is the shift of the peak en- 
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Figure 4. The observed ES{E) spectra of the prompt emission 
of GRB 980425 (top panel) as observed by BATSE and of GRB 
031203 (bottom panel) as observed by INTEGRAL, expressed in 
fluence units. The data shown for GRB 980425 have been an- 
alyzed by us using a power law with an exponential cut-off to 
fit the time integrated spectrum over the ~ 30 sec duration of 
the burst. For GRB 031203, we report the spectrum analyzed by 
us, using version 5.1 of the OSA analysis software and the cor- 
responding response matrices. This spectrum corresponds to the 
first 20 sec of the emission as the following 20 to 40 sec do not 
contribute significantly. Solid lines corresponds to the flux passing 
through a screen of material, while dashed lines are the intrinsic 
flux (i.e. before passing through the screen). The parameters of 
the model are listed in Table 1. 



ergy -Bpeak to higher values, as shown in Fig.|21 This is due 
to the declining (with energy) Klein-Nishina cross section 
(7KN- For instance, the transmitted spectrum in the case of 
tt = 6 has -Epeak ~ 2 MeV, while the intrinsic spectrum has 
Speak = 511 keV. 

Note also that the bolometric transmitted energy, for 
this particular case of an intrinsic -Bpeak = 511 keV, de- 
creases much less than the factor ea:p(— rx) expected in the 
Thomson regime (a factor ~10 for tt ~ 5, instead of the 
factor 150 we would have in the Thomson regime). 



3.2 Application to GRB 980425 and GRB 031203 

In Fig.^Jwe show the observed spectrum of the two sources 
(time integrated over the duration of the bursts), the as- 
sumed intrinsic spectrum, and the transmitted one. 

We reanalyzed the INTEGRAL data obtained with the 
IBIS/ISGRI instrument using version 5.1 of the OSA anal- 
ysis software and the corresponding response matrices. The 
spectrum corresponding to the first 20 seconds of the burst 
(start time 2003-12-03 22:01:26) can be fit by a power law 
with photon index 1.7 and fiuence 1.4e-6 erg cm~^. 

The found slope is consistent with the slope found by 
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Figure 5. The correlation between -Epeak a^nd the isotropic 
equivalent energy radiated during the prompt emission Ei^^ (blue 
and green symbols, the "Amati" relation, Amati et al. 2002)) and 
the coUimation corrected energy E^ = (1 — cos (?j )-Eiso, where 6j 
is the semiaperture angle of the jet, assumed conical (red crosses, 
the "Ghirlanda" relation, Ghirlanda et al. 2004, but here shown 
in the case of a circumburst density profile n oc appropri- 
ate for a stellar wind made by the progenitor, see Nava et al. 
2006). Blue points are the 19 GRBs with a firm estimate of the 
jet opening angle; green points are those bursts (adapted from 
Amati et al. 2006) without a firm estimate of 9j . GRB 980425 
and GRB 031203 are outliers (open circles) with respect to the 
"Amati" relation (long— dashed line), but they can be made con- 
sistent with this correlation dotted line and green filled circles) by 
assuming that the radiation we see is the transmitted radiation 
crossing a scattering screen as illustrated in Fig. |3] Without the 
scattering screen, these two bursts would have appeared as inter- 
mediate between X-ray rich GRB and X-ray flashes. Moreover, 
they can be made consistent with the Ghirlanda relation ( solid 
line, in the case of a "wind" density profile) if their jet opening 
angle is 9° and 11°, as labelled. GRB 060218 would be consistent 
with the Ghirlanda relation for 8j ~ 21°. Alternatively we also 
show the effect of spectral evolution (Sec. 4) which can make the 
two outliers consistent with the "Amati" relation (short-dashed 
line and blue-filled squares) if we only detected a short portion of 
their much longer emission which rapidly evolved from the hard 
7— ray band to the soft X-ray band. In this case, if we require that 
these two bursts are consistent with the "Ghirlanda" relation we 
derive angles of 34° and 13° for GRB 980425 and GRB 031203, 
respectively. 



Sazonov et al. (2004), but smaller by a factor 1.4 in normal- 
isation. 

For GRB 980425 we have analysed the BATSE data (to 
show the data points) using a cut-off power law model. The 
best fit parameters are a = 0.45 ± 0.22 and -Epeak ~ 138 
keV (x^ ~ 109/97 dof), giving a flux of ~ 10~^ erg/cm'^ sec 
(40-700 keV), consistent with that reported in Plan et al. 
2000. 

The same model (power law -I- exponential cut) was 
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GRB 


7-hp 




F, 


peak 

keV 


^obs 
peak 

keV 


ft- 


f 
J 


980425 
031203 


7.5 
6.5 


0.7 
0.85 


83.3 
200 


25 
30 


138 
292 


11° 
9° 


577 
161 



Table 1. Input parameters of the scattering model for the two 
outUers. The parameter / is the ratio of the intrinsic to trans- 
mitted total energy. The jet opening angle 9j is derived requiring 
that the the intrinsic iJpeak a-^d Ej fit the "Ghirlanda" relation 
derived under the wind circum-burst density hypothesis (Nava et 
al. 2006). 

assumed for GRB 031203, for which we plot the spectrum 
analyzed by us in Fig. 2] but with frequencies in the source 
rest frame. For GRB 031203 a soft component has been in- 
ferred through the observed transient rings produced by the 
scattering of the burst radiation by dust layers in our own 
Galaxy (Vaughan et al. 2004, Tiengo & Mereghetti 2006 - 
see also Vaughan et al. 2006 for another case of dust scat- 
tering halo associated with GRB 050724). This component, 
that will be discussed in the next section, is here neglected. 

It can be seen in Fig. |1] that the model can reproduce 
the observed data quite satisfactorily. In Tab. 1 we report 
the assumed parameters (tt, a, -Be) for GRB 980425 and 
GRB 031203, and in Fig.|^we show where these two bursts 
"move" in the i?pcak--Eiso plane (filled circles) when the inci- 
dent spectrum (i.e. that impinging on the scattering screen) 
is reconstructed with the assumed parameters. In Fig.|K|we 
also show what jet opening angle Oj they should have to also 
fit the "Ghirlanda" relation between Speak and the coUima- 
tion corrected energy E-y = (1 — cos &j)ii'iso, as derived under 
the hypothesis of a wind density environment (Nava et al. 
2006). This angle turns out to be ~ 9° for GRB 031203 and 
~ 11° for GRB 980425. These values of 9j lie in the large- 
angle tail of the distribution of known jet opening angles 
(for 19 GRBs - see Fig. 6 in Nava et al 2006). In fact, the 
low peak energies of the intrinsic spectra (i.e. before pierc- 
ing through the scattering screen) places both GRB 980425 
and GRB 031203 in the XRF region of the Speak--B plane 
of Fig. 13 where the "Amati" and "Ghirlanda" correlation 
converge. Bursts in this region have, on average, jet open- 
ing angles larger than those placed at the opposite extreme 
(such as e.g. GRB 990123, the most energetic burst in the 
Epea]i-EiBo and in the Ep^ak-Ej planes). 

Although this "screen model" might have interesting 
feature, it suffers from some problems, related to the pres- 
ence of an optically thick screen itself. The first problem is 
why this material is present in some bursts but not in others. 
An additional problem concerns its location: if it is at some 
distance from the bursts, its material would not be ionized 
(and the dust in it not destroyed) by the prompt radiation, 
making the optical afterglow not detectable. If it is instead 
located close to the burst, one should see some signature of 
its presence when the fireball and/or the supernova ejecta 
collide with it. 



4 SPECTRAL EVOLUTION 

The scenarios presented in the previous sections were aimed 
to explain the prompt emission as observed by BATSE and 



BeppoSAX for GRB 980425, and by INTEGRAL for GRB 
031203, assuming that the received prompt emission was 
mainly in the 20-300 keV range. This may not be true, how- 
ever, for the following reasons: 

• The echo-rings detected by XMM-Newton a few hours 
after the trigger of GRB 031203 can be produced only if 
the incident flux is at soft X-ray energies and large. Re- 
constructing the flux, fluence and spectrum of the emission 
which was scattered by the dust layers of our Galaxy is not 
straightforward, but the existing estimates point towards a 
total fluence at least as large as the fluence detected by IN- 
TEGRAL, at energies around 1 keV (Watson et al. 2006; 
Tiengo & Mereghetti 2006). In this respect, note that the 
re-analysis of the INTEGRAL data performed by us, with 
improved detector response matrices yields a total fluence a 
factor 1.4 less than what quoted by Sazonov et al. (2004), 
while the spectral shape is the same. 

• One important issue is if this soft component belongs 
to the prompt emission or if it is instead the beginning of 
the afterglow. The large fluence, however, suggests that this 
is part of the prompt, since we never observed X-ray after- 
glows as bright as the prompt in any other GRB. 

• Another issue concerns the duration of the soft com- 
ponent. A lower limit can be roughly estimated based on 
the the lack of detection with IBIS/ISGRI in the images ob- 
tained after (and before) the GRB prompt emission. Such 
images are limited to the energy range above ~17 keV. 
Therefore one has to invoke either that the soft X-ray emis- 
sion desumed from the dust scattering rings had a spectrum 
with a very sharp cut-off before the IBIS/ISGRI energy 
range or that the soft X-ray emission lasted more than sev- 
eral hundred seconds, in order to give the required fluence 
without violating the flux upper limits. Note that a dura- 
tion of the order of a few thousand seconds would still be 
compatible with the observed width of the expanding rings 
measured with XMM-Newton. 

• For GRB 980425 we do not have direct evidences of 
the presence of a (possibly long-lasting) soft component, 
but we do have some indication of spectral softening from 
the analysis of the WFC (Wide Field Camera) and GRBM 
monitor on board BeppoSAX, as presented by Frontera et al. 
(2000), see the spectra reported in Fig. |^ The total fluence 
in the WFC [2-28 keV] is of the same order than the fluence 
in GRBM (which agrees with BATSE). 

• The idea of an hard-to-soft spectral evolution, with a 
very long soft emission, received a strong support from the 
recently Swift GRB 060218 (Cusumano et al. 2006). For this 
burst the narrow field instruments of Swift could follow the 
emission for most of the prompt phase, which lasted at least 
^3000 seconds (Barthelmy et al. 2006). This burst showed 
a hard-to-soft spectral evolution, and its time integrated 
spectrum has a peak energy around 5 keV and a ~ 6 x 10'*^ 
erg isotropic energy which makes this burst fully consistent 
with the Amati relation (Campana et al. 2006; see Fig. 1^. 
The total energetics of GRB 060218 is very similar to the one 
of GRB 031203, and both are associated with a supernova of 
similar power (e.g. Cobb et al. 2006). This burst prompted 
us to explore in some detail the idea of spectral evolution 
for the two outliers. 

• Finally we point out that, recently, a tight correlation 
was discovered (Firmani et al. 2006) , involving only observ- 
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ables belonging to the prompt emission: the peak luminosity, 
Speak and the "high signal" timescale To. 45 (it is the same 
timescale used for the characterization of variability, see Re- 
ichart et al. 2001). GRB 980425 and GRB 031203 are outliers 
with respect to this correlation. However GRB 031203 would 
become consistent with this newly found correlation if -Bpcak 
were in the ~keV range, as suggested by the echo-rings. 

In this section we see if it is possible to construct a 
simple model for the spectral evolution of GRB 031203 and 
GRB 980425 that is at the same time consistent with the 
existing observations of these two GRBs and also predicting, 
for their time integrated spectra, values of -Epcak and iJiso in 
agreement with the Amati relation. We are inspired to do 
such a modelling by the behaviour of GRB 060218, and we 
will then first construct our toy model in order to reproduce 
the time-dependent behaviour of its spectra. For this burst 
we have much more information than for the others, and 
the modelling therefore merely implies to characterise the 
spectral index, the cut-off energy and the normalisation of 
the spectrum with smooth functions of time, to reconstruct 
both the time resolved spectra and the light curves in the 
XRT [0.3-10 keV] and BAT [15-150 keV] energy bands. 

For GRB 031203 the time dependent information is al- 
most not existent: we assume that the soft X-ray spectrum 
and fluence, around 1 keV, is is the one inferred (albeit in 
an approximate way) , by the echo-rings observed by XMM- 
Newton. 

For GRB 980425 we have even less information, and 
the modelling is done only to show that there is no evidence 
against a hard-to-soft evolution also for this burst. 

We are aware that i) this model does not pretend to 
be unique, since it can be one of many spectral evolution 
models; ii) it is completely phenomenological, and it has no 
theoretical support, iii) for the sake of simplicity, we will 
assume that the spectral indices and the high energy cut 
off of the spectrum, as well as its normalisation, are smooth 
functions of time (i.e. we will assume no flares), and finally 
iv) even if we will minimise the number of free parameters, 
in the end we need many. 

These can be taken as strong limitations if we were re- 
ally seeking "the" model for interpreting the properties of 
spectral evolution, but this is not our aim. Our goal is in- 
stead to demonstrate that, even with a naive and simple 
approach, it is indeed possible to fulfil the two requirements 
of being in agreement with all observations and nevertheless 
making these specific bursts to obey the Amati relation. 

4.1 The case of GRB 060218 

Fig. |S] (upper panel) shows the time integrated spectrum 
of GRB 060218. We re-analyzed the XRT data of the first 
~2600 sec (in wt mode) starting ~ 159 sec after the BAT 
trigger (03:34:30 UT - Cusumano et al. 2006). For the data 
extraction we used the standard xrtpipeline (v. 0.9. 9) and for 
the spectral analysis we used the v. 007 detector response 
matrix (for wt mode limited to grades 0-2). Campana et al. 
2006, through a joint XRT-BAT fit, found that the peak 
energy -Epcak ~ 4.5 keV lies in the XRT energy range, and 
the spectrum can be fitted by the sum of a black-body (BB) 
and a cut-off power law (GPL). We show in Fig.|S|the X-ray 
data as re-analysed by us, together with the best fit model 



(dashed lines). The duration of this burst was so long that 
the BAT "burst" mode data contained only the first ~ 350 
sec of its emission. The study of the following portion of 
the burst requires the analysis of the "survey" mode data 
which are more difficult to handle, we believe that the BAT 
team can analyse these data more appropriately than us. 
However, the BAT team already gives the total BAT fiuence 
for this burst (Campana et al. 2006). We have then fitted 
the XRT time integrated spectrum with the same model (i.e. 
BB-I-CPL and same A''h = 6 x lO^'^ cm~^ - see Campana 
et al. 2006), using the BAT fluence (Campana, priv. comm.; 
Dai, Zhang & Liang 2006) as a constrain. In this way we 
have obtained the same -Epcak and 0.3-10 keV fluence of 
Campana et al. (2006). 

Through this "BAT-fiuence-constrained" fit procedure 
we analyzed the XRT data in the same time intervals re- 
ported in Campana et al. (2006). Fig. |S| shows three spec- 
tra: the top one is the spectrum integrated over the total 
observed duration of the bursts, while the other two spectra 
correspond to the [159-309 s] and [2456-2748 s] time bins. 
These are the first and the last spectra that can be obtained 
with both XRT and BAT data, according to the time bin- 
ning chosen by Campana et al. (2006), which can be also 
seen in the right bottom panel of Fig. |S| For the [159-309 
s] time-bin the BAT "burst" mode data were available and 
we can then show the corresponding XRT-I-BAT data. For 
the [2456-2748 s] time-bin we only show the XRT data, fit- 
ted with the same model and procedure (i.e. forcing the fit 
to have the same BAT fluence reported in Campana et al. 
2006) as for the time-integrated spectrum. 

Note also that for all the three spectra (first and last 
time bin and time integrated) we show the BB component 
resulting from the fit (dotted lines). This component con- 
tributes at least 20% to the total spectral fiuence, and it has 
been interpreted as due to the jet shock breakout (Cam- 
pana et al. 2006). For this, we will not consider it in our 
spectral-evolution model. 

In the right bottom panel of Fig. |S| we show the light 
curves in the BAT and XRT energy ranges and also the 
evolution of Ec (i.e. Ec = i?pcak/(l — ct), where a is the 
energy spectral index of the CPL model). The values of Ec 
have been derived with the same procedure discussed above. 

To model the time evolution of this burst we need the 
normalisation K and the spectral index a of the power law, 
and the cut-off energy E^, as a function of time. For the 
latter we have simply interpolated the derived values with 
a smooth function (the combination of a very weakly rising 
and a decaying power law of time: solid line in the right 
bottom panel of Fig. |SJ, while the assumed behaviour of 
K and a are shown in the left bottom panel of the same 
figure. We have chosen functions which are either constants 
or power laws of time. 

Having reconstructed in this way the entire time evolu- 
tion of this burst, we could then sum-up the instantaneous 
spectra to give the spectra integrated in specific time bins 
and the spectrum integrated over the entire burst duration. 
Note that the time integrated spectra are not necessarily de- 
scribed by a BB+CPL model, since the spectral parameters 
evolve. This can be the reason why, quite often, one obtains 
a better fit with a time resolved spectrum rather than with 
a time integrated one. 

The agreement with the observation is excellent, but 
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Figure 6. Top panel: Spectra of GRB 060218 for different time-bins: i) entire duration (top); ii) [159-309 s] (rising spectrum with 
also BAT data) iii) [2456-2748 s] (soft spectrum). We plot ES{E) vs E, S{E) being the fluence. Dotted lines indicate the blackbody 
component, not considered for the spectral evolution, and long— dashed lines represents the best fit obtained from the analysis of the data. 
Continuous lines show the results of our proposed modelling. Left bottom panel: assumed behaviour of the normalisation K and energy 
spectral index a. Right bottom panel: light curves in the BAT (15—150 keV) and XRT (0.3-10 keV) range, and evolution of E^. The 
flux in the 0.3—10 keV is the (de-absorbed) flux of the cut-off power law component only: we have subtracted the blackbody component 
from the total flux. Continuous lines are the results of our modelling. 



expected, since what we have done is merely to interpo- 
late with smooth functions of time the real evolution of the 
spectral parameters. The overprediction of the flux at high 
energies in the time integrated spectrum (solid-red line in 
top panel) is due to the fact that the smooth function inter- 



polating Ec lies above the data points in the second part of 
the time evolution (see the right bottom panel of Fig. |SJ . 
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Figure 7. Top panel: the spectral energy distribution of GRB 031203 is compared with the result of our model. Here we show the 
INTEGRAL data as analysed by us and the low energy emission as inferred by the light scattering echo calculated by Watson et al. 
(2006). The bow-tie corresponds to the uncertainty in the amount of the scattering material. In addition to those, there are additional 
uncertainties connected to which cross section is used, as discussed in Tiengo & Mereghetti (2006). The lines corresponds to the emission 
time integrated in different energy intervals, as labelled. The bottom left panel shows how the normalisation Kit), the high energy 
spectral index a and the cut-off energy Ec change in time. The right bottom panel compares the observed and the calculated light curves 
in the three labelled energy bands. 



4.2 GRB 031203 

This burst was discovered by INTEGRAL, and the spectrum 
in the 20-400 keV band is shown in Fig. |7| In this figure we 
show the spectrum as reanalysed by us (see. Section 3.2). 



The soft X-ray spectrum is inferred by the light scatter- 
ing echo discovered by XMM-Newton (Watson et al. 2004) 
and modelled by Watson et al. (2006) to get the intrin- 
sic spectrum illuminating the galactic dust clouds. Note 
that Tiengo & Mereghetti (2006) re-derived the illuminating 
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spectrum and fluence by adopting a different cross section 
for the scattering process between X-ray photons and dust, 
obtaining a fluence somewhat smaller than Watson et al. 
(2006), but with a similar spectrum. Therefore both groups 
agrees that the soft (less than a few keV) X-ray emission lies 
much above the extrapolation from the high energy spec- 
trum, and that the slope of the spectrum below a few keV 
is softer than the 20-200 keV slope. 

The lightcurves (in counts per second) in tree differ- 
ent energy bands have been extracted by us using the most 
recent version of the response matrix. 

The modelling of the spectral evolution of this burst is 
aimed to reproduce: i) the time integrated spectrum as seen 
by INTEGRAL in the 20-200 keV band; ii) the upper hmit 
in the medium energy band [17-40 keV]; iii) the soft X-ray 
emission as inferred by the light scattering echo and iv) the 
light curves (in counts per second) in three different energy 
bands of INTEGRAL [20-50 keV; 50-100 keV and 100-200 
keV]. 

4-2.1 Modelling the spectral evolution 

We model the spectral evolution of the prompt emission as- 
suming that the spectrum, at any time, is a broken power 
law ending with an exponential cut. We assume a broken 
power law (instead of a single one) to limit the total ener- 
getics to a finite value. For simplicity, we assume that the 
energy spectral index a of the low energy branch is equal 
to 1/2 the value of the energy spectral index of the high 
energy part (we do this just to limit the number of free pa- 
rameters). The two power laws connect at the break energy 
Eh- We assume that the cut-off energy Ec, the break energy 
Eh, the power law index a and the normalisation can evolve 
in time. We set 

^(t)/2 



F{E,t) = K(t)\ — 



F{E,t) 



rp s -ait) 



exp 



E-Eh 
' E4t) 



E<Eh 



E>Eh 



(8) 



where the normalisation K is 

here 6i and &2 are indices describing the evolution of the 
normalisation of the spectrum, and tm is related to the time 
imax of the maximum of K{t) by 

/bi \ 1/(61+62) 

tmax = tm(-r] (10) 



The spectral index a and the cut-off energy, Ec, are assumed 
to evolve as power laws of time (i.e. a oc t"^ , Ec oc t^'^, but we 
allow for the possibility that they remain initially constant 
(for a time equal or shorter than imax). All these choices are 
arbitrary, the only guiding line is simplicity. 

The constant K is chosen to match the observed spec- 
tral shape of the fluence the burst obtained by the high 
energy data. 

To reconstruct the light curves in one band AE = E2~ 
El keV in counts we multiply F{E, t) by the normalised 
effective area A{E) of the INTEGRAL instruments to get 




10 100 
E [keV] 

Figure 8. Calculated spectra of GRB 031203 at different times, 
as labelled (number are seconds). Also shown is the time inte- 
grated spectrum (time interval 0—3000 seconds). 



C{AE,t) = 



E2 



El 



F{E,t) 
E 



A{E)dE 



(11) 



Then we re-normalise all the light curves obtained in differ- 
ent energy bands by the same normalisation factor. 

4-2.2 Results for GRB 031203 

The left bottom panel of Fig.|7]shows the behaviour of K{t), 
a{t) and the cutoff energy Edt) which we have assumed for 
the spectral evolution of GRB 031203. 

The top panel of Fig. Q shows the SED resulting from 
our modelling, while the right bottom panel of the same 
flgure reports the light curves of the model (solid lines) in 
different bands and compares those with the data. The as- 
sumed overall duration of the prompt emission is 3000 sec- 
onds: of those, only the first ~ 50 seconds have been de- 
tectable by INTEGRAL in the 20-200 keV energy band. 

We can see that the assumed hard-to-soft spectral evo- 
lution is in agreement with all the information we have about 
the prompt emission of this burst, and predicts a peak en- 
ergy of the overall spectrum at -Bpcak ~ 5 keV. Also the total 
energetics is enhanced if we include the soft X-ray prompt 
emission. This is shown in Fig. |51 where one can see how 
GRB 031203 "moves" in the Biso-^Epcak plane. With the 
derived values of these parameters GRB 031203 obeys the 
Amati relation. 

4.3 GRB 980425 

The prompt emission of this burst was detected both by 
BATSE and by BeppoSAX. The latter detected it both with 
the Gamma Ray Burst Monitor (GRBM) and the Wide 
Field Camera (WFC) (e.g. Plan et al. 1999). While in Fig. 
13] we show the BATSE spectrum (as re-analysed by us), in 
the top panel of Fig.|5|we show the BeppoSAX data (Fron- 
tera et al. 2000). This is because we are here interested in 
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Figure 9. Top panel: the spectral energy distribution of GRB 980425 is compared with the result of our model. Here we show the 
BeppoSAX data, as originally reported by Frontera et al. (2000), in three different time— intervals, and the spectrum time-integrated 
over the duration of the burst as seen by BeppoSAX. The bottom left panel shows how the normalisation Kit), the high energy spectral 
index a and the cut-off energy Ec change in time. The right bottom panel compares the observed and the calculated light curves in the 
three labelled energy bands. The observed light curves are from BATSE. 



the medium energy X-ray band more than the high energy 
one. The time integrated fluences of BATSE and the GRBM 
instrument are consistent. We also show the spectra time- 
integrated in 3 different time intervals, using the data pub- 
lished in Frontera et al. (2000). 

In the right bottom panel of Fig. |2l we show the hght 
curve of the count rate as detected by BATSE, in three 



different energy bands [25-50 keV; 50-100 keV and 100-200 
keV]. 

We have then applied our toy model using the same 
approach described in Sec. 4.2.1. The evolution of K(t) a(t) 
and i5c(t) is shown in the left bottom panel of Fig. 1^1 

The resulting spectra are shown as solid lines in the 
top panel of Fig. |^ They have been time integrated in the 
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Observed 

-^pcak 
Eisa 


GRB 980425 
138 

i.DZX iU 


GRB 031203 
291 

i.zX iU 


keV 

erg 


Scattering 








-^pcak 


25 


30 


keV 




y.oo X lu 


l.yo X lU 


erg 


Spectr. Evo. 

-^pcak 
Elso 


0.6 

2.54x1048 


5 

1.6x10^0 


keV 

erg 



Table 2. The value of 

^peak and -Eiso that are observed and 
that are predicted by the scattering and spectral evolution models 
discussed in this paper. These are the values used in Fig|5] Note 
that the Siao (observed) for GRB 031203 takes into account the 
small reduction due to the re-analysis of the INTEGRAL data. 

same time intervals as the data. We also show the spectrum 
corresponding to the complete time evolution, assumed to 
last for 3000 seconds. We have chosen this duration just 
because it is equal to the one assumed for GRB 031203 and 
observed for GRB 060218. The calculated light curves are 
overlapped to the BATSE 64ms observed light curves in the 
right bottom panel of Fig. |5| 

We can see that, again, the calculated spectra and light 
curves are consistent with what observed. For this burst, 
unfortunately, we do not have any supplementary informa- 
tion forcing us to believe that the duration of the prompt 
emission is long (i.e. the WFC was not sensitive enough to 
follow the burst for a long time, and there was no informa- 
tion about late emission as was the case of the "echo rings" 
for GRB 031203). 

On the other hand, our point is the following: if GRB 
980425 is similar to GRB 031203 (and both are similar to 
GRB 060218), then it is possible that its peak energy -Epoak 
is much smaller than what derived by the high energy spec- 
tra only. In the case we have just shown -Epcak ^ 0.6 keV, 
once the spectrum is integrated over the entire (assumed) 
duration. Were this the case, then GRB 980425 would not 
be an outlier any longer, as can be seen in Fig. |3 where we 
show the new location of this burst. 



5 DISCUSSION 

In this work we have explored three different scenarios to 
see if GRB 980425 and GRB 031203 are really outliers with 
respect to the Amati relation (and therefore also outliers to 
the Ghirlanda relation) or if they only appear as such. 

The first model (off axis viewing angle) can bring the 
two bursts on the Amati relation only if the energetics of 
these bursts (as would be seen by an on-axis observer) is 
very large. Such bursts are not expected to exist at low red- 
shifts, and then this argues against this model as an expla- 
nation for the two bursts being only "apparent" outliers. 

Then we explored the idea that these bursts have "nor- 
mal" energetics and values of iJpoak, but both these quanti- 
ties are modified by a scattering cloud located in the vicinity 
of the bursts. It must be located in the vicinity of the burst 
because in this case the prompt emission can completely 
ionise the material (Lazzati & Perna 2002), making the soft 



X-ray afterglow insensitive to absorption. If very close to the 
burst, the scattering material would be accelerated to large 
velocities. In the comoving frame of the moving material, 
incoming photons would be seen redshifted, and a larger 
fraction of photons would scatter in the Thomson regime. 
The transmitted spectrum would then be similar (albeit not 
identical) to the one obtained in the case of a not moving 
material with a larger tt. The complete and exact analysis 
implies to assume the location and width of the scattering 
screen, which would determine its acceleration, and this lies 
outside the aim of this paper. 

In this model the intrinsic (i.e. before scattering) ener- 
getics is greater than the observed one, but not by a large 
factor. Even with a substantial value of the scattering op- 
tical depth (i.e. rx ~a few), the reduction in flux of the 
transmitted radiation is of the order of one-a few hundreds. 
The scattering process is more effective on low energy X-ray 
photons, implying that the transmitted spectrum is "bluer" 
than the intrinsic one. Therefore the intrinsic -Epcak can be 
smaller than observed. The two things together (i.e. intrin- 
sically, the two bursts are more energetic and much redder) 
make it possible for them to lie on the Amati relation. The 
increase in energetics is not so large to face the problem of 
having too energetic bursts at small redshifts. 

The problem, for this idea, is to explain why these 
bursts are surrounded by dense and thick material (for scat- 
tering), while other bursts are not. The partial answer we 
can give is to consider one obvious selection effect: we see 
distant bursts only if they are surrounded by material with 
a small value of tt, since otherwise they go under the de- 
tection threshold. Consider that even values of tt ~ 1 have 
almost no effect on the observed spectrum, since in this case 
more than half of the photons pass through the scattering 
cloud undisturbed (see Fig. 3). Even the fast variability, if 
present, would not be smeared out, since the received flux is 
dominated by the transmitted photons: the scattered ones, 
having a longer path to travel, arrive later and more diluted 
in time (i.e. with a much reduced flux). 

On the other hand, this model poses several other prob- 
lems. The most severe of them is that we expect some sig- 
nature of the presence of such a thick screen in the vicinity 
of the burst, produced by the collision of the bursts and 
supernova ejecta with the screen itself. 

The third model we have proposed is the most promis- 
ing, especially after GRB 060218, which we used as a guide 
to model the two outliers. In this framework the derived in- 
trinsic powers and energetics are the least demanding, as can 
be seen in Fig. |3 and in Table 2, where we list the values of 
-Eiso and iSpcak for the two bursts. In the framework of this 
model the presence of both a hard and a soft component (in 
GRB 031203) is the result of the spectrum evolving in time, 
and is not due to two separated components (i.e. two emis- 
sion mechanisms). Note that the real bolometric and time 
integrated energetics (i.e. Eiso) is larger than what derived 
by considering the high energy spectrum only, but not by a 
large factor. 

The main result of our study is that GRB 031203 and 
GRB 980425 are likely to be apparent outliers: their intrin- 
sic properties are instead consistent with the Amati rela- 
tion (and therefore it is possible that they obey even the 
Ghirlanda relation, and they would do so with jet opening 
angles not particularly extreme, see Fig. |K1 although larger 
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than the average values found for "normal" GRBs). This 
possibility is strengthened by the fact that GRB 060218, 
energetically a twin of GRB 031203, lies on the Amati rela- 
tion independently of our modelling. GRB 980425 and GRB 
031203 could also become consistent with the newly found 
correlation between the peak luminosity of the prompt, 
Speak and the "high signal" timescale To. 45 (Firmani et al. 
2006), if they indeed have iJpcak around ~1 keV. 

On the other hand these two bursts, together with GRB 
060218, besides being underluminous, have a duration much 
longer than the other bursts, including the X-ray rich and 
the X-ray flash population of bursts (see e.g. Lamb, Don- 
aghy & Graziani 2005) . The inevitable conclusion is that the 
Amati relation (and possibly the Ghirlanda one) is more ro- 
bust than previously thought, being obeyed by bursts that 
are underluminous and with a duration longer than the aver- 
age one and that (with respect to the Amati relation) GRB 
980425 and GRB 031203 are not representative of a different 
population of objects with respect to "classical" GRBs. 
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